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Abstract 

Complete one-loop results for the decay widths of neutral Higgs bosons (h a ) into lighter 
neutral Higgs bosons (hf,, h c ) are presented for the MSSM with complex parameters. 
The results are obtained in the Feynman-diagrammatic approach, taking into account 
the full dependence on the spectrum of supersymmetric particles and all complex phases 
of the supersymmetric parameters. The genuine triple-Higgs vertex contributions are 
supplemented with two-loop propagator- type corrections, yielding the currently most 
precise prediction for this class of processes. The genuine vertex corrections turn out to 
be very important, yielding a large increase of the decay width compared to a prediction 
based on the tree-level vertex. The new results are used to analyse the impact of the 
experimental limits from the LEP Higgs searches on the parameter space with a very 
light MSSM Higgs boson. It is found that a significant part of the parameter space 
of the CPX benchmark scenario exists where channels involving the decay hi — > h\h\ 
have the highest search sensitivity, and the existence of an unexcluded region with 
Mfn ~ 45 GeV is confirmed. 
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1 Introduction 



Higgs self-couplings, i.e. triple-Higgs couplings, h a hbh c , and quartic Higgs couplings, h a hbh c hd, 
are a crucial element of electroweak symmetry breaking via the Higgs mechanism. While the 
prospects for a direct experimental determination of the quartic Higgs coupling at present 
and future colliders are small (see e.g. Ref. [1]), probing the triple-Higgs coupling will be 
one of the prime goals in the experimental programme for testing the Higgs mechanism. 
This coupling can be accessed via a precision measurement of the Higgs production pro- 
cess e + e _ — > Zh a h a at the ILC [2] or CLIC [3], and via Higgs cascade decays of the form 
h a — * h b h c . While Higgs cascade decays are obviously impossible in the Standard Model 
(SM), they can play an important role in models with extended Higgs sectors. 

Besides the interest in Higgs cascade decays as a means to directly probe Higgs self- 
couplings, a precise prediction for decays of this kind is also important for phenomenological 
reasons. Where kinematically possible the Higgs cascade decay modes can even be dominant 
and affect Higgs phenomenology very significantly. 

A well-known example of an extended Higgs sector is the Minimal Supersymmetric Stan- 
dard Model (MSSM), whose Higgs sector comprises two scalar doublets that accommodate 
five physical Higgs bosons. In lowest order these are the light and heavy CP-even h and H, 
the CP-odd A, and the charged Higgs bosons H ± . Higher-order contributions yield large 
corrections to the masses and couplings, and also induce CP-violation leading to mixing 
between h, H and A in the case of general complex SUSY-breaking parameters. The corre- 
sponding mass eigenstates are denoted as h\, h 2 , hs. If the mixing between the three neutral 
mass eigenstates is such that the coupling of the lightest Higgs boson to gauge bosons is 
significantly suppressed, this state can be very light without being in conflict with the exclu- 
sion bounds from the LEP Higgs searches [5]. In this case the second-lightest Higgs boson, 
h 2 , may predominantly decay into a pair of light Higgs bosons, h 2 — > hih\. 

The results of the Higgs searches at LEP [4, 5] have been interpreted in certain MSSM 
benchmark scenarios [6,7]. In the CPX scenario [7], which involves large complex phases 
of the trilinear coupling, A t , and the gluino mass parameter, M3, the decay h 2 — > h\h\ is 
important in a significant part of the parameter space. This decay mode leads to a quite 
complicated final state, corresponding to a 6-jet topology if h 2 is produced in association 
with an hadronically decaying Z boson or in association with h\. In the analysis of the LEP 
Higgs searches within the CPX scenario an unexcluded region for light h\ and relatively 
small values of tan (3 (the ratio of the vacuum expectation values of the two Higgs doublets) 
remained, so that no firm lower bound on the mass of the lightest Higgs boson of the MSSM 
could be set. The unexcluded parameter region with a very light Higgs boson will also be 
difficult to cover with the Higgs searches at the LHC [1,8,9] (see also Ref. [10] for a recent 



In order to reliably determine which parameter regions of the MSSM with a very light 

1 It should be noted that Higgs cascade decays are possible in the MSSM also in the limit where complex 
phases are neglected. The decay h — > AA occurs in a small parameter region with very light Ma [11], leading 
to small unexcluded parameter regions in the Ma -tan /3 plane from LEP Higgs searches [5] (especially in 
the "no- mixing" scenario [6]). For large values of Ma also the decay H — > hh can occur. Higgs cascade 
decays into very light Higgs bosons can also be important in a considerable part of the parameter space of 
extensions of the MSSM, see e.g. Ref. [12] for a discussion within the NMSSM. 
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Higgs boson are unexcluded by the Higgs searches so far and which regions will be accessible 
by Higgs searches in the future, precise predictions for the Higgs cascade decays h a — > h h c 
in the MSSM with complex parameters (cMSSM) are indispensable. For propagator-type 
corrections the evaluations are quite advanced, and results incorporating the dominant 
one- and two-loop contributions have been obtained within the Feynman-diagrammatic 
(FD) approach [13-15] and the renormalisation-group (RG) improved effective potential 
approach [16]. The public codes FeynHiggs [14,15,17,18], based on the FD approach, and 
CPsuperH [19], based on the RG improved effective potential approach, are available. For 
the genuine h a hbh c vertex contributions, on the other hand, so far only effective coupling 
approximations have been available in the cMSSM. 

In this paper we obtain complete one-loop results within the FD approach for the decays 
h a — > hbh c taking into account the full dependence on all complex phases of the supersym- 
metric parameters. The mixing of the three neutral Higgs bosons among themselves and 
with the Z boson and the unphysical Goldstone-boson degree of freedom are consistently 
taken into account. We furthermore obtain complete one-loop results for the decays of neu- 
tral Higgs bosons into SM fermions in the cMSSM, h a — > //. The new one-loop results 
are combined with all existing higher-order corrections in the FD approach, yielding in this 
way the currently most precise predictions for the class of processes h a — > ht,h c . The re- 
sults presented in this paper will be included in the code FeynHiggs [14,15,17,18]. We 
find that the genuine vertex corrections are very important for a reliable prediction of the 
Higgs cascade decays. The genuine vertex corrections lead to a drastic change compared 
to a prediction taking into account propagator-type corrections only. We compare our full 
result with various approximations and investigate the dependence on the complex phases. 
As an application of our improved theoretical predictions, we analyse to what extent the 
previously unexcluded parameter region with a rather light Higgs boson is compatible with 
the limits on topological cross sections obtained from the LEP Higgs searches. 

2 The MSSM with complex parameters: notations and 
conventions 

The MSSM, in its most general form, contains CP-violating phases at tree level in the Higgs, 
slepton, squark, chargino/neutralino and gluino sectors. The gauge-boson, lepton and quark 
sectors do not contain extra phases (we assume in this paper a unit CKM matrix). In our 
calculation the full dependence on all complex phases of the supersymmetric parameters is 
taken into account. In the following we briefly specify the notations and conventions used 
in this paper and define the parameters that are relevant for the discussion of the numerical 
results. 

2.1 The squark sector 

The bilinear terms of the squarks in the Langrangian give rise to the mass matrix 
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where 
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A q — //{cot /3, tan/3}, 
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and the trilinear couplings A q and the Higgs-mixing parameter \x can be complex. This mass 
matrix needs to be diagonalised to get the tree- level physical states qi,q 2 , 
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where in our conventions Cq is real and Sq is complex. 

2.2 The neutral Higgs sector 
2.2.1 Tree level 

We write the two Higgs doublets as 
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where v\ and v 2 are the vacuum expectation values, and tan/3 = v 2 /v±. The conventions 
used here are the same as in Ref. [14]. The MSSM Higgs sector is CP-conserving at lowest 
order. The tree-level neutral mass eigenstates h, H, A and the unphysical Goldstone-boson 
degree of freedom G are related to the CP-even neutral fields 0i, 02 and the CP-odd neutral 
fields Xi? X2 through a unitary matrix, 



(5) 



In the renormalisation prescription that we have adopted, the parameter tan/3 receives a 
counterterm, while the mixing angle (3 Q remains unrenormalised. After the renormalisation 
has been carried out, one can set (3d — (3. 

2.2.2 Higgs mass matrix in higher orders 

To find the loop-corrected neutral Higgs masses, the poles of the 3x3 propagator matrix 
A(jo 2 ) in the (h, H, A) basis need to be found (mixing with the Goldstone boson G and the 
Z boson can be neglected in the propagator matrix since the corresponding contributions 
are of sub- leading two-loop order, see the discussion in Ref. [14]). Determining the poles of 
the propagator matrix is equivalent to finding the three solutions to 
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The propagator matrix is related to the 3x3 matrix of the irreducible renormalised 2-point 
vertex- functions F 2 (p 2 ) through 

[-Afr 2 )]- 1 = f 2 (p 2 ) = i [p 2 l-M(p 2 )] , (7) 

where 

/ml-Z hh (p 2 ) -t hH {p 2 ) -t hA {P 2 ) \ 

M(p 2 ) = -t hH {p 2 ) ml-t HH {p 2 ) -t HA {p 2 ) . (8) 
V -t hA (f) -£ha{p 2 ) m\-t AA {p 2 )j 

Here m^, m#, m A are the lowest-order mass eigenvalues, and the S^- = h,H,A) are 
the renormalised self-energies. In general, the three solutions M- 2 ha (with h a = hi, h 2 , h 3 ) are 
complex. They can be written as 

M 2 ha = Ml-iM ha W ha , (9) 

where Mh a is the (loop-corrected) mass of the respective Higgs boson, Wh a is its width, and 
by convention 

M hl < M h2 < M h3 . (10) 
We calculated the self-energies using an expansion about M 2 a , 

%k{M 2 ha ) » t jk (M 2 J + zlm [M 2 ha ] t' jk (Ml), (11) 

with j, k = h,H,A (as a check of our procedure, each time Yijk(M.\ a ) was calculated, the 
next term in the expansion was also explicitly calculated, to ensure it is negligible). To find 
each solution to Eq. ([H]), an iterative procedure was used. 

3 Complete one-loop calculation of Higgs-boson cas- 
cade decays and decays into SM fermions 

We calculated the full 1PI (one-particle irreducible) one- loop vertex contributions to the pro- 
cesses h a — > hbh c , taking into account all sectors of the MSSM and the complete dependence 
on the complex phases of the supersymmetric parameters. Examples of generic diagrams for 
one of the contributing topologies are shown in Fig. [TJ 

In order to obtain precise predictions for the branching ratios BR(/i a — > hbh c ) it is im- 
portant to calculate also the decay widths of the Higgs bosons into SM fermions, h a — > f f, 
at the same level of accuracy, since the decay modes into bb and r + r~ are dominant over 
large parts of the cMSSM parameter space. We therefore derived also complete one-loop 
results for the processes h a — > // (including SM-type QED and, where appropriate, QCD 
corrections) for arbitrary values of all complex phases of the supersymmetric parameters. 
The partial widths for the other Higgs-boson decay modes have been taken from the program 
FeynHiggs [14,15,17,18]. 

Our calculations have been carried out in the FD approach, making use of the pro- 
grams FeynArts [20] and FormCalc [21]. Concerning the renormalisation, we use the same 
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Figure 1: Examples of generic diagrams (showing only one of the topologies) contributing to 
the processes hi — ► hjhf., where hi, hj, = h, H, A. Furthermore, / are SM fermions, / are 
their superparters, x°, \ are neutralinos and charginos, V are vector bosons, H denote the 
neutral and charged Higgs bosons and the Goldstone bosons, u are Faddeev-Popov ghost 
fields. 



transformations and renormalisation conditions as in Ref . [14] . We parameterise the electric 
charge in the lowest-order decay amplitudes in terms of oc(Mz), corresponding to the charge 
counterterm 
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Here II 7 (0) is the photon vacuum polarisation, S T denotes the transverse part of the self- 

(5) 

energy, and Aa = Aa^ ad + Aai cpt is the shift in the fine-structure constant arising from large 
logarithms of light fermions. The other parameter renormalisations are listed in Ref. [14]. 
The fermion fields in the processes h a — > ff are renormalised on-shell. For the renormalisa- 
tion of the Higgs fields it is convenient to use a DR scheme as in Ref. [14], while the correct 
on-shell properties of the S-matrix elements involving external Higgs fields are ensured by 
the inclusion of finite wave-function normalisation factors. 
The wave-function normalisation factors are obtained from 
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where the A hh (p 2 ), A hH (p 2 ) etc. are the elements of the 3x3 Higgs propagator matrix (see 
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Ref. [14] for expressions in terms of renormalised self-energies), and h a , hj,, h c are some per- 
mutation of hi, h 2 , h 3 . For the evaluation of self-energies at complex momenta the expansion 
in Eq. fTTTj) is employed. 

The wave-function normalisation factors can be expressed in terms of a (non-unitary) 
matrix Z [14], 



(16) 




where 



(17) 



Here f\ a is a one-particle irreducible n-point vertex-function (including loop corrections) 
which involves a single external Higgs boson h a . 
The matrix Z fulfills the conditions 

lim - \ (Z-f 2 -Z T ) 

V^M\ a P 2 -M 2 ha \ Jhh 



lim - - \ (z-f 2 -Z T ) = 1 (19) 
v^M\ b P 2 -M 2 hb \ Jhh 

lim - - \ (z-t 2 -Z T \ = 1, (20) 
P ^Ml c p 2 -M 2 hc \ Jaa 

where r 2 has been introduced in Eq. ([7j). We choose in the following h a — h±, hi, — h 2 
and h c = h%. It should be noted that this choice is purely a convention. Other choices 
would give the same physical results. This fact can most easily be seen if the wave-function 
normalisation factors are defined as in Eqs. (fT3l) - (|T5l) . i.e., with the Higgs-boson self-energies 
evaluated at the complex pole. In this way the evaluation of the masses and the wave-function 
normalisation factors is treated on the same footing. The definition of the wave-function 
normalisation factors adopted in Eqs. (TT3l) - (115j) differs slightly from the one in Refs. [14,22], 
where the wave-function normalisation factors were defined at the real part of the complex 
pole (furthermore, in Refs. [14, 22] the real parts of the diagonal wave-function normalisation 
factors Zh, Zh, Za were taken). The two definitions of the wave-function normalisation 
factors differ by contributions from imaginary parts that are formally of sub-leading two- 
loop order. It turns out that the inclusion of the imaginary parts improves the numerical 
stability of the results in certain parameter regions, while otherwise these contributions are 
completely negligible (see also Ref. [23]). 

In a complete one-loop calculation of Higgs decay processes also mixing contributions 
between the states h, H, A and the Goldstone and Z bosons have to be taken into account 
(we denote these reducible contributions as r^™"). We treat these mixing contributions 
strictly at one-loop level in order to ensure the cancellation of unphysical poles (while our 
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prescription for the wave function normalisation factors described above automatically incor- 
porates leading reducible higher-order contributions). Since the mixing self-energy involving 
G and Z is already a one-loop contribution, in a strict one-loop expansion the r^™ x are 
evaluated at the (unrotated) tree- level masses m 2 h .. Accordingly, our results can be written 
as (summation over repeated indices is understood) 



■pfull 

-pfull _ 
1 haff 



plPI 

hihjh k 



r S/«)+ r r 



m h^ m h 



m 



hk 



(21) 
(22) 



where as before hi, hj, hk = h, H, A. The numerical impact of the mixing contributions with 
the Goldstone and Z bosons on the results presented in this paper turned out to be small. 

In our numerical analysis below we will compare our full result with the contribution 
obtained from the t, t sector in the approximation where the gauge couplings are neglected 
and the diagrams are evaluated at zero external momenta. We refer to this approximation 
as the "Yukawa Approximation", which is expected to yield the leading one- loop contribu- 
tion if tan j3 is not too large. In this approximation, the counterterm contributions in the 
renormalised vertex all vanish. 



4 Combination with higher-order results 

We have combined our new one-loop result with the most up-to-date higher-order propagator- 
type contributions in the FD approach. This has been done by supplementing the one-loop 
self-energies in Eqs. flH]) and (|T3|) - (|T5l) with the two-loop self-energies obtained from the 
program FeynHiggs [14,15,17,18]. In FeynHiggs the 0(a>ta s ) corrections are incorporated 
including the full phase dependence at the two-loop level, while other two-loop contributions 
so far are only known in the limit of vanishing complex phases [24] . In our numerical analysis 
below we will restrict to those higher-order contributions for which the phase dependence at 
the two-loop is explicitly known, i.e. we do not include the residual contributions for which 
a result exists only in the MSSM with real parameters. The calculation of the decay width 
r(^a - * bb) furthermore contains resummed SUSY-QCD corrections, including the full phase 
dependence. 

By supplementing our complete one-loop results for the processes h a — > hbh c and h a — > // 
with the state-of-the-art higher-order propagator-type corrections we obtain the currently 
most precise predictions for the h a — > hf,h c decay widths and branching ratios in the MSSM 
with complex parameters. It should be noted that in the special case where the complex 
phases are put to zero our results also provide improved predictions for the Higgs cascade 
decays occurring in the CP-conserving case, i.e. h — > AA and H — * hh. The numerical 
impact of the latter will be discussed elsewhere. 

In our numerical analysis we also investigated the impact of using loop-corrected Higgs 
masses and couplings (rather than the tree-level values) within loop diagrams, which is for- 
mally a higher-order effect. The numerical impact on our results turned out to be negligible. 
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5 Implementation of exclusion bounds from the LEP 
Higgs searches 

A precise prediction for the process hi — > h\h\ is particularly important for investigating 
the exclusion bounds from the Higgs searches at LEP [4, 5] within the MSSM with complex 
parameters. The LEP analysis in the CPX benchmark scenario [7] resulted in an unexcluded 
parameter region for relatively small tan (3 where the lightest neutral Higgs boson is very light 
but has strongly suppressed couplings to gauge bosons. The second-lightest Higgs boson, on 
the other hand, may be within the kinematic reach of LEP in this region but has a large 
branching ratio into a pair of lightest Higgs bosons, i.e. hi —>■ h\h\. 

In our numerical analysis below we analyse the impact of our new result on the LEP 
exclusion bounds in the cMSSM. This is done by comparing the cMSSM predictions with 
the topological cross section limits given in Refs. [4,5,25]. The LEP limits on the various 
cross sections [4, 5] have been implemented into the code HiggsBounds [26] (for applications 
of preliminary versions of HiggsBounds, see Ref. [27]). In order to obtain a correct statistical 
interpretation of the overall exclusion limit at the 95% C.L., the first step is to determine, 
for each parameter point, which one of the various channels has the highest statistical sensi- 
tivity for setting an exclusion limit [25]. We then compare the theoretical prediction for this 
particular channel with the topological cross section limit determined at LEP for this chan- 
nel. We neglect, in this analysis, the theoretical uncertainties from unknown higher-order 
corrections — see Ref. [26] for a discussion of this issue. The predictions of the topological 
cross sections have been obtained using the wave-function normalisation factors defined in 
Sect. El 

While the topological cross section limits are very convenient for testing a wide class of 
models and are not restricted to specific parameter values, it should be kept in mind that 
the dedicated analyses carried out in Ref. [5] for certain benchmark scenarios have a higher 
statistical sensitivity than the limits obtained from the topological cross sections. This is 
in particular the case in parameter regions where two (or more) channels have a similar 
statistical sensitivity, since the method based on the topological cross section limits allows 
one only to use one channel at a time. 



6 Numerical results 

In our numerical analysis we use the parameter values of the CPX benchmark scenario [7] , 
adapted to the latest experimental central value of the top-quark mass [28] and using an 
on-shell value for the absolute value of the trilinear couplings A t and that is somewhat 
shifted compared to the DR value specified in Ref. [7] (see Ref. [29] for a discussion in the 
MSSM with real parameters). Specifically, if not indicated differently we use the following 
parameters (the lowest-order Higgs-sector parameters tan/3 and M^± are varied; in our 
analysis we include tan (3 values up to 40 and M H ± values up to 1000 GeV.) 

Msusy = 500 GeV, \A t \ = \A b \ = 900 GeV, fi = 2000 GeV, m~ g = \M 3 \ = 1000 GeV, 
mt = 170.9 GeV, M 2 = 200 GeV, (23) 
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Figure 2: Full result for T(h 2 — > hih{) compared with various approximations, see text, (a) 
r(/i2 h\hi) vs. M/ ll for tan/5 = 6 (M#± is varied; all other parameters are set to the CPX 
values given in Eqs. ([23]), (El)), (b) T(h 2 -> Mi) as function of tan/3 for M hl = 30 GeV 
(M H ± is adjusted to ensure = 30 GeV). 



and the complex phases of the trilinear couplings A t , and the gluino mass parameter M3 
are set to 

7T 7T 7T . . 

<PA t = -, <PA b ^9= 2" ( 2A ) 

In Eq. (|23|) Msusy denotes the diagonal soft SUSY-breaking parameters in the sfermion mass 
matrices that are chosen to be equal to each other, Msusy = Ml = Mg R , see Eq. (TjQ). 



6.1 Results for the hi —> h\hi decay width 

Fig. [2] shows the relative effect of different contributions to the h 2 — > h\h\ decay width in 
the area of the cMSSM parameter space which is particularly relevant to an investigation of 
the unexcluded regions in the LEP Higgs searches [5]. In Fig. [2] (a) the full result for the 
decay width and various approximations are plotted against M hl while keeping tan/5 = 6. 
All the results plotted include the higher-order corrected wave-function normalisation factors 
as described in Sects. [3] and HI They differ only in the genuine contributions to the h 2 h\h\ 
vertex. One can see that the full result (denoted as 'Full') differs drastically from the case 
where only wave-function normalisation factors but no genuine one-loop vertex contributions 
are taken into account ('Tree'). The inclusion of the genuine vertex corrections that have 
been evaluated in this paper can increase the decay width by more than a factor of six in 
this example (for values of sufficiently below the kinematic limit of M/ ll = 0.5Mh 2 where 
the decay width goes to zero). The Yukawa approximation agrees much better with the full 
result, giving rise to deviations of up to ~ 30%. Using the full contribution from the t, t, b, b 
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Figure 3: r(/i 2 — » hihi) as function of tan/3 for M/ ll = 30 GeV and different values of tfA t 
(Mjj± is adjusted to ensure M hl = 30 GeV). 

sector ('t, st, b, sb') and from all three generations of fermions and sfermions ('f, sf ') yields 
a prediction that deviates from the full result by up to 20%. 

In Fig. [2] (b), the decay width is plotted against tan/3 whilst adjusting M H ± such that 
M/ ll = 30 GeV. The meaning of the various lines is the same as in Fig. [2] (a). The result 
where the genuine vertex corrections are neglected ('Tree') shows a completely different 
qualitative behaviour and differs drastically from the full result. The pronounced dependence 
on tan/3, giving rise in particular to a region where T(/i 2 — » hihi) ~ for tan/3 ~ 4.3, is due 
to the fact that the Yukawa vertex corrections (which are dominant) to the matrix element 
change sign when tan (5 is varied. As the decay width depends on the absolute value squared 
of the matrix element, the region where the Yukawa vertex corrections are close to zero 
corresponds to a minimum of the decay width. The deviations between the full result and 
the contribution from only the fermion and sfermion sector reach about 15% in this example. 

In Fig. [3]the full result for T(/i 2 — ¥ h\hi) is shown as a function of tan (3 for M/ ll = 30 GeV 
(adjusting M H ± accordingly) and different values of <pA t - The dependence on <pA t is ver y 
pronounced, leading in particular to a relative shift of the curves in tan j3. As a consequence, 
comparing the results for r(/i 2 — * h\h\) for different values of (fA t at fixed values of tan /3 can 
yield dramatic effects. Thus, a thorough treatment of the phase dependence is indispensable 
for a meaningful theoretical prediction of Higgs cascade decays in this parameter region. 

6.2 Analysis of exclusions bounds from the LEP Higgs searches 

In Figs. H] and [5] the M/jj-tan/3 parameter space of the CPX scenario is analysed. Fig. H] 
shows the branching ratio of the Higgs cascade decay, BR(/i 2 — > h\hi). It can be seen that, 
over a large part of the parameter space where the decay /i 2 — > h\h\ is kinematically possible, 
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Figure 4: The branching ratio BK(h 2 — > hihi) in the M^-tan/3 plane of the CPX scenario. 

it is actually the dominant decay channel. The branching ratio is particularly large for low 
and moderate values of tan (3. In the region where tan (3 ~ 4-5 the Yukawa contribution to 
the matrix element for the h 2 —>■ h-Jii decay changes sign and causes a sharp drop in the 
h 2 — > h\h\ branching ratio, as was already observed in Fig. [2j A similar behaviour occurs 
also in the region of tan f3 ~ 35. 

Plot (a) of Fig. [5] indicates which channel has the highest statistical sensitivity and there- 
fore which channel will be used to set an exclusion limit in different regions of the parameter 
space. As explained in Sect. [5], this information is needed for an interpretation of the topolog- 
ical cross section limits obtained at LEP [4,5] as 95% C.L. excluded regions in the M^-tan/3 
parameter space. One can see in the figure that the channels e + e~ — > h 2 Z — > h\h\Z — > bbbbZ 
and e + e~ — > h 2 hi — > hihihi — > bbbbbb have the highest search sensitivity in a region with 
small Mh l and moderate values of tan j3, 5 ^ tan (3 % 9. For small M/ ll and somewhat 
higher tan/? the channel e + e~ — > (hihi) — > (bbbb) has the highest search sensitivity. It 
should be noted that all channels involving the decay of the h 2 boson in the region of small 
M/ ll are strongly influenced by the T(h 2 — > hihi) decay width, either directly in the case of 
the channels involving the Higgs cascade decay, or indirectly through the branching ratio of 
the h 2 . The parameter region where T(h 2 — * h\hi) is important coincides with the region 
of the CPX scenario that could not be excluded at the 95% C.L. in the analysis of the four 
LEP collaborations [5]. 
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Figure 5: Coverage of the LEP Higgs searches in the M^-tan/3 plane of the CPX scenario. 
Plot (a) shows the channels that are predicted to have the highest statistical sensitivity for 
setting an exclusion limit. The colour codings in are: red (■) = (h\Z) — > (bbZ), blue (■) = 
(h 2 Z) -> (JbbZ), white (□) = (h 2 Z) ->_{h x h x Z) -> (bbbbZ), cyan (■) = (hJn) -> (6656), 
yellow ( ) = (h 2 hi) — > (hih±hi) — > (bbbbbb), black (■) = other channels. Plot (b) shows the 
parameter regions excluded at the 95% C.L. by the topological cross section limits obtained 
at LEP [4,5]. The colour codings are: green (dark grey) = LEP excluded, white = LEP 
allowed. 



In plot (b) of Fig. [5] we have combined our new theoretical predictions (containing the 
complete one-loop result for the genuine vertex corrections and higher-order corrections, 
as described above) with the topological cross section limits obtained at LEP. We find an 
unexcluded region at Mh 1 ~ 45 GeV and moderate tan f3 where channels involving the decay 
h 2 — > hihi play an important role. Thus, our analysis, based on the most up-to-date theory 
prediction for the h 2 —* h\h\ channel, confirms the 'hole' in the LEP coverage observed in 
Ref. [5] (see in particular Fig. 19 of Ref. [5]). 

It should be noted, on the other hand, that the results shown in plot (b) of Fig. [5] differ 
in several respects from the results presented in Ref. [5]. As discussed above, near to borders 
between areas where different search topologies are predicted to have the highest exclusion 
power our analysis has less statistical sensitivity than the benchmark studies of Ref. [5]. A 
further difference is the input value of the top-quark mass. While we are using the latest 
experimental central value of m t = 170.9 GeV [28], most of the analysis of Ref. [5] was done 
for m t = 174.3 GeV. We have explicitly checked that (as expected) the unexcluded region 
is significantly increased if we use m t = 174.3 GeV instead. Concerning differences in the 
theoretical predictions, the analysis of Ref. [5] was based on the two codes FeynHiggs2. 0, an 
earlier version of the FeynHiggs program [18], and CPH [7], a predecessor of the program 
CPsuperH [19]. For each scan point, in Ref. [5] the results from CPH and FeynHiggs2.0 (for 
the decay width T(h 2 — > h-Jix) the CPH formula was used in both codes) were compared 
with each other, and the result yielding the more conservative exclusion bound was retained. 
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Figure 6: Impact of the genuine vertex corrections on the 'hole' in the LEP coverage for 
M hl w 45 GeV and tan (3 ~ 7. Plot (a) shows the full result (detailed view of plot (b) of 
Fig. [5]). Plot (b) shows the result for the case where only contributions from SM fermions 
and their superpartners are taken into account in the genuine vertex corrections. Plot (c) 
shows the result where the Yukawa Approximation has been used for the genuine vertex 
corrections. The colour codings are: green (dark grey) = LEP excluded, white = LEP 
allowed. 



Our theoretical prediction for T{]i2 — * hxh\) based on a complete diagrammatic one-loop 
result of the genuine vertex contribution goes significantly beyond the effective coupling 
approximation used in Ref. [5]. Furthermore, the 0(a t a s ) propagator-type corrections ob- 
tained in the Feynman-diagrammatic approach for arbitrary complex phases [15] were not 
yet available when the analysis of Ref. [5] was carried out. 

In Fig. [6]we focus on the uncovered parameter region at M/ ll ph 45 GeV and tan /3 ~ 7 and 
compare our complete one-loop result for the genuine triple-Higgs vertex corrections with two 
approximations. While plot (a) shows the full result (i.e., it is a detailed view of plot (b) of 
Fig. ED, in plot (b) the genuine vertex corrections are approximated by the contributions from 
fermions and sfermions only, and in plot (c) the Yukawa Approximation has been used for 
the genuine vertex corrections. In all three plots the wave function normalisation factors and 
all other decay widths are the same (calculated as described above). The differences between 
the three plots in Fig. [6] are therefore entirely caused by the genuine vertex corrections to 
the Higgs cascade decays. While all three graphs show an unexcluded region at M/ ll m 
45 GeV and tan ph 7, the shape of this region changes quite significantly. In particular, 
the analysis based on the full result gives rise to a considerably larger unexcluded region 
around ~ 45 GeV and tan/5 ~ 7 compared to the Yukawa Approximation (while the 
Yukawa Approximation gives rise to a larger unexcluded region for smaller M^J . As expected 
from Fig. [21 the inclusion of all SM fermion and sfermion loop contributions yields a better 
approximation of the full result. 



7 Conclusions 

We have obtained, within the MSSM with complex parameters, complete one-loop results for 
the classes of decay processes in which a heavier neutral Higgs boson decays into two lighter 
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ones and for those in which a neutral Higgs boson decays into a pair of SM fermions. Our 
results take into account all sectors of the MSSM and incorporate the full dependence on 
all complex phases of the supersymmetric parameters and the external momenta. The new 
one-loop results have been supplemented with state-of-the-art propagator-type corrections 
containing higher-order contributions obtained in the FD approach (taken from the program 
FeynHiggs), yielding in this way the currently most precise predictions for these processes 
within the MSSM with complex parameters. 

We find that the genuine vertex contributions to the triple-Higgs vertex are numerically 
very important. Their inclusion changes the predictions for the decay widths very drastically 
compared to the case of an effective coupling approximation based only on the wave function 
normalisation factors of the external Higgs bosons. Including genuine vertex corrections in a 
simple Yukawa Approximation yields a prediction for the decay width that is much closer to 
the full result, but we still find deviations of up to 30% in the numerical examples investigated 
in this paper (using the CPX scenario). Even contributions beyond the fermion/sfermion 
sector can have a sizable impact on the decay widths. We have furthermore found that the 
dependence of the decay width on the complex phase tfA t is very pronounced, emphasizing 
the need for a thorough treatment of the effects of complex phases. 

Based on our improved theoretical predictions we have analysed the impact of the limits 
on topological cross sections obtained from the LEP Higgs searches on the parameter space 
with a very light Higgs boson within the MSSM with complex parameters. We find that, over 
a large part of the parameter space where the decay h 2 — > h\h\ is kinematically possible, it is 
the dominant decay channel. The corresponding search channels e + e~ — > h 2 Z — > h\h\Z — > 
bbbbZ and e + e~ — > h 2 hi — > h-Jiihi — > bbbbbb have the highest sensitivity for setting an 
exclusion limit in a region with small M/ ll and moderate values of tan (3. We find that a 
parameter region with ~ 45 GeV and tan/3 »s 7 remains unexcluded by the limits 
on topological cross sections obtained from the LEP Higgs searches, confirming the results 
of the four LEP collaborations achieved in a dedicated analysis of the CPX benchmark 
scenario. A precise theory prediction for the h 2 — > h\hi channel is crucial for mapping out 
the unexcluded parameter region (it should be noted in this context that in the parameter 
region of the CPX scenario also formally sub-leading two-loop corrections can have a sizable 
numerical impact). We find that the shape of the unexcluded region is significantly modified 
if the full result for the vertex corrections is replaced by approximations. 

The results presented in this paper will be included in the public code FeynHiggs. It would 
be interesting to compare our results with the other public code for evaluating Higgs masses 
and decay widths in the MSSM with complex parameters, CPsuperH. This comparison is 
affected, however, not only by the genuine triple-Higgs vertex corrections that are the main 
focus of the present paper, but also by differences in the propagator-type corrections used 
in the two codes. Furthermore, a meaningful comparison between FeynHiggs and CPsuperH 
requires a translation between the on-shell input parameters used in FeynHiggs and the DR 
input parameters used in CPsuperH. Such a detailed comparison is beyond the scope of the 
present paper. We will address this issue in a forthcoming publication. 
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